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We investigate a silicon single-electron transistor (SET) in a metal-oxide-semiconductor (MOS) 
structure by applying a magnetic field perpendicular to the sample surface. The quantum dot is 
defined electrostatically in a point contact channel and by the potential barriers from negatively 
charged interface traps. The magnetic field dependence of the excitation spectrum is primarily driven 
by the Zeeman effect. In the two-electron singlet-triplet (ST) transition, electron-electron Coulomb 
interaction plays a significant role. The evolution of Coulomb blockade peaks with magnetic field B is 
also owing to the Zeeman splitting with no obvious orbital effect up to 9 T. The filling pattern shows 
an alternate spin-up-spin-down sequence. The amplitude spectroscopy allows for the observation 
of the spin blockade effect, where the two-electron system forms a singlet state at low fields, and 
the spin polarized injection from the lead reduces the tunneling conductance by a factor of 8. At a 
higher magnetic field, due to the ST transition, the spin blockade effect is lifted and the conductance 
is fully recovered. 

PACS numbers: 73.23.Hk, 72.25.-b, 73.40.Qv 



I. INTRODUCTION 

As a promising candidate of spin qubits for quan- 
tum computing, semiconductor quantum dots have at- 
tracted tremendous research effort .'i' Recent progress in 
GaAs lateral quantum dots has demonstrated all pos- 
sible single-qubit operation^ and the square-root-of- 
swap operation on two qubits,'^ which form a universal 
set for quantum computing.^ Compared to GaAs, lat- 
eral quantum dots in silicon are expected to have a 
spin coherence time orders of magnitudes longer, be- 
cause ^^Si has no nuclear spin, and there is no hy- 
perfine interaction between electron spins and nuclear 
spins.l^^ There has been steady progress in the devel- 
opment of silicon-based single electron transistors us- 
ing silicon-on-insulator (SOI) wafers,^ Si/SiGe quantum 
well structures® and metal-oxide-semiconductor (MOS) 
structures,'^ though it remains a considerable challenge 
due to material properties.!^ 

In this paper, we present our recent progress on an 
enhancement-mode MOS single electron transistor (SET) 
in pure silicon.^ The unique structure of the MOS-SET 
allows us to verify the formation of quantum dots and 
better understand the device, which is largely neglected 
in some early works. When properly biased, the inves- 
tigated SET has a quantum dot in an electrostatically 
defined point contact channel. In addition to single elec- 
tron tunneling behavior, we have also observed its mag- 
netic field dependence. For all of the data presented here, 
the field is applied perpendicular to the sample surface. 
The observed magnetic field dependence of the excitation 
spectruirP is found to be primarily driven by the Zeeman 
splitting. Furthermore, the spectrum enables us to di- 
rectly observe the singlet-triplet (ST) transition, where 
electron-electron Coulomb interaction plays a significant 
role. The evolution of Coulomb blockade peaks with the 



magnetic field is also measured. The data strongly sug- 
gest that the ground state energies also shift with the 
applied magnetic field owing to the Zeeman effect. Up to 
9 T, there is no obvious orbital effect. The evolution of 
peak amplitudes illustrates the spin blockade effect due 
to magnetic field induced spatial separation of spin-up 
and spin-down states. At a higher magnetic field, the 
spin blockade effect is lifted when the SET undergoes 
the ST transition. 



II. SAMPLES AND TRANSPORT 
CHARACTERISTICS 

Figures [ija) andjljb) show the schematic top view and 
the cross-sectional view of the device. All devices are fab- 
ricated on N-type, high purity silicon (100) wafers with a 
resistivity of 3 ^ 5 kfi-m. By using photolithography and 
ion implantation, the heavily phosphorous-doped regions 
are defined as the source (1) and the drain (2). A 27 nm 
thick thermal oxide is grown by dry oxidation. It also 
has a bilayer-gated configuration. Six side gates (labeled 
as A-F in Fig.jlja)) are located above the thermal oxide 
and buried in the second dielectric layer, Si02, which is 
400 nm in thickness and grown by high-density plasma- 
enhanced chemical vapor deposition (HDPECVD). The 
top gate (Gl), which is above the HDPECVD oxide, 
laterally overlaps with the ohmic source and drain re- 
gions. The device is annealed in forming gas at 420 °C 
for 30 minutes to reduce the interface states. The fab- 
rication process and the basic operating principle have 
been discussed in more detail elsewhere.'^ These devices 
are characterized in a dilution refrigerator with ^10 mK 
base temperature. The electron temperature is about 
400 mK by fitting the Coulomb blockade peaks with the 
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FIG. 1: (Color online) (a) Schematic top view and (b) cross 
sectional view of a MOS-SET. In (a), the shaded regions, 
the solid lines, and the rectangular area, depict the heav- 
ily phosphorus-doped ohmic source (1) and drain (2) leads, 
the 6 side gates (A-F), and the top gate (Gl) respectively. 
It also includes an equivalent circuit with two small quan- 
tum dots, due to local potential fluctuation, and one large 
quantum dot at the center, which is electrostatically defined 
by the 6 side gates, (c) Scanning electron micrograph of 
70 nm-wide side gates before the top Si02 layer is grown 
by high-density plasma-enhanced chemical vapor deposition 
(HDPECVD). The gap, d, between two neighboring side gates 
is ~160 nm. The circle depicts the location of the quantum 
dot of interest as discussed in the text, (d) Transfer character- 
istics of single quantum dots. The gate biases are: Vci = 16 
V, Vb = Vf = 1 V. In the first trace (depicted as Va,b), 
the side gate voltages Va,b {~ Va = Vb) are swept, when 
Vc = Vd = V. In the second trace (depicted as Vc,d), 
the side gate voltages Vc,d (= Vc = Vd) are swept with 
Va = V^B = V. 

equation jliSl 

G/G,nax = cosh-\ea{Vg - Vga)/2kBT), (1) 

where G is the source-drain conductance, Gmax is the 
maximum conductance at the peak, VgQ is the gate volt- 
age at resonance, ks is the Boltzmann constant, and T 
is the electron temperature. The ratio of the gate ca- 
pacitance to the total capacitance, a = Gg/CY,, is ex- 
tracted from the slopes of nearby diamonds in the stabil- 
ity chart. The source-drain conductance is measured by 
the standard ac lock-in technique using a 37 Hz, 0.1 mV 
excitation voltage. The positively biased top gate (Gl) 
induces two-dimensional (2D) electrons at the silicon and 
silicon thermal oxide interface. The peak electron mobil- 
ity is about 5200 cni'^/Vs at 400 mK in a test hall bar 
device, which has the same structure, but without the 6 
side gates. When ^ Vb = Vc = Yd = Ve = Vf ^ ^ 
V, the SET device turns on at Vqi > 9 V. (Transport 
data not shown here.) The negatively biased side gates 
A and B can deplete electrons below and define a point 
contact channel. Due to the local potential fluctuation 
near the Si/Si02 interface, there is a quantum dot formed 



in the channel, as illustrated in Fig. [ija). Side gates C 
and D, and separately, gates E and F, can also be used 
to define a quantum dot through the same mechanism. 
When the gates are properly biased, the device displays 
single quantum dot SET characteristics. As shown in 
Fig. [T]^d), two typical data illustrate the transfer behav- 
ior. In the first trace (depicted as V^,b), Vqi = 16 V, 
Ve = Vf = 1 V, = Vi) = V, and side gates A 
and B {Va = Vb = Va^b) are biased from to —2.5 V. 
As Va.b is swept, the drain conductance oscillates; when 
Va^b < —2.2 V, the conductance is less than 1.3 x 10~^ 
(e^/h), the noise fioor of our measurement system. In the 
second trace (depicted as Vc,d), the gate voltage Vc,d 
(= Vc = Vd) is swept, while Va — Vb — V, and the 
SET conductance also oscillates as Vc^d < —1.3 V and 
becomes diminished at —2.5 V. These pictures are well 
supported by the following data. 

To verify the formation of the quantum dots and to 
better understand the operation of the device, we have 
investigated the dependence of the source-drain conduc- 
tance 6*21 on Va,b and Vc.d- The observed transfer 
characteristic shown in Figj2la), a portion of which is 
enlarged and shown in Fig. pp)) , demonstrates a double 
quantum dot behavior. When Va,b and Vc,d are both 
less than about —1.2 V, the device is in the weakly cou- 
pled double quantum dot regime. In this weakly coupled 
regime, two triple points will merge into ones, i.e. the 
high conductance spots as clearly shown in Fig. Jgtb), 
and are located at the vertices of parallelograms .'^How- 
ever, when —1.2 V < Vc^d < 1 V, sweeping Va^b re- 
produces the single electron tunneling features similar to 
that shown in Fig. [ijd). The same is observed when 
sweeping Vc,d while keeping —1.2 V < Va,b < 1 V. 
The nearly horizontal and vertical lines in Fig. [2]ja) in- 
dicate that these two dots are spatially separated. One 
dot is strongly capacitively coupled to side gates A and 
B, but only weakly coupled to gates C and D. The other, 
however, is strongly coupled to gates C and D, and only 
weakly coupled to gates A and B. That is, one dot is 
physically near side gates A and B, whereas the other 
dot is located close to side gates C and D. Furthermore, 
data shown in Fig. [2|c) suggest that the quantum dot 
is located in the narrow point contact channel. Keeping 
the other two point contacts fully conductive {Vci = 16 
V, Ve = Vp = I V, and = Vb = V) and sweeping 
the voltage of side gates C and D in the range of interest, 
the source-drain conductance shows clear single electron 
tunneling characteristics. Because the capacitances be- 
tween the dot and, separately, side gates C and D are 
the same, the peak positions in Fig.|2jc) display diagonal 
dependence.!^ To be more specific, the quantum dot is 
physically located in the point contact channel, at equal 
distances to side gates C and D. The discontinuity of the 
diagonal lines is caused by a mere one electron charge 
change at a nearby charge trap. Our experiment shows 
that this nearby charge trap has a negligible effect on our 
results in this paper. A similar characteristic is also ob- 
served for the dot defined by side gates A and B. (Data 



3 




1HIHII0IIN1-2 1 M-^ 




V{x) V{x) 

FIG. 2: (Color online) Characteristics and equivalent circuits 
of the device at Vai = 16 V, Vb = Vp = 1 V. (a) Source-drain 
conductance 6*21 as a function of side gate voltages Va,b and 
Vc,D- (b) Enlarged section of (a) shows the underlining di- 
agonal lines from the center large quantum dot. (c) Source- 
drain conductance G21 as a function of Vc and Vd , with 
Va = Vb ~ V- Diagonal lines suggest that side gate C and 
side gate D are coupled to the small quantum dot with the 
same capacitance, as illustrated in (e). (d) Equivalent circuit 
and schematic potential profile under multiple quantum dots 
bias conditions, (e) Equivalent circuit and schematic poten- 
tial profile under single quantum dot bias conditions. 



not shown here.) These observations are consistent with 
the picture that there are two weakly-coupled quantum 
dots, located in the point contact channels defined by 
side gates A, B, and C, D, respectively. Residue inter- 
face charges are known to form potential fluctuation at 
the Si/Si02 interface.'^^It has been reported earlier that 
a point contact channel in silicon, due to these i nterface 
charges, can show Coulomb blockade oscillations .'^^Ulllil 
In addition to the two quantum dots discussed above, 
side gates A, B, C, D, E, and F, can also be biased to 
form a (large) quantum dot. A closer look at Fig. l2Fa) 



indeed reveals the underlining diagonal lines from this 
dot, as shown in Fig. [2]jb). 

In summary, under proper bias condition, Va.b < — 1-2 
V and Vc,D < —1-2 V, the equivalent circuit consists of 
an SET at the center (with a large quantum dot), in 
series with two others (both with small quantum dots), 
as illustrated in Fig. [2|d); on the other hand, a pair of 
side gates can define a single SET, when the other side 
gates are positively biased, as shown in Fig.|2]^e). These 
pictures do capture the major characteristics of the de- 
vice. However, there is not always a quantum dot in a 
point contact. When the gap between two neighboring 
side gates is reduced from ~ 160 nm to ~ 90 nm, the 
point contacts can smoothly turn off the device without 
Coulomb oscillations in some of our samples. In these 
devices, the electrostatically defined dot is the only fea- 
ture, but they currently only work in the many-electron 
regime. With continuous down-scaling and improvement, 
it is quite possible to reach the few-electron regime, even 
realize single electron confinement. 



III. SINGLE QUANTUM DOT AND 
MAGNETIC FIELD SPECTROSCOPY 

In the following, for the purpose of investigating single 
spin in silicon, we focus on the single dot SET confined 
in the channel defined by side gates C and D. Fig. [sj^a) 
shows the stability chart of the SET, where the source- 
drain differential conductance G21 is measured against 
the source-drain dc bias, V21, and the side gate voltages 
Vc,D, while Ve = Vf ^ 1 Y and Va = Vb = V. The 
observation of Coulomb diamonds confirms the forma- 
tion of a single quantum dot in the point contact channel. 
Within the orthodox theory, the quantum dot is modeled 
by a disc with a diameter d, and the total capacitance 
Cs is Aed, where e (= 11.9 in silicon) is the dielectric 
constant. From the diamond shown in Fig. |3ja), we can 
directly measure the half height (V21 — e/Cs), and the 
obtained charging energy Ec {— e^/Cs) is about 6 meV. 
So Cs is approximately 27 aF, which suggests a disc di- 
ameter of about 60 nm. This estimated disc diameter is 
of the same order of magnitude as the averaged spacing 
between the interface-trapped charges in our device; we 
have separately measured by the conductance method 
using large area devices the interface quality, and ob- 
tained, at room temperature, an interface trap density 
of ~ 6 X 10^" eV^^cm^^, corresponding to a spacing be- 
tween charges to be ~ 40 nm. The origin of the quan- 
tum dot of interest is probably coming from negatively 
charged electron traps, located at the Si/Si02 interface. 
Single dopants such as phosphorous atoms had been pro- 
posed to be behind transport studieJ^^l in silicon SETs, 
but our observed addition energy and the excitation en- 
ergy are both much less than what one would expect 
using the single dopants picture. In addition, as will be 
discussed below, we have observed a spin filling pattern 
different from that of the single phosphorous dopants pic- 
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FIG. 3: (Color online) (a) Stability chart taken with Vai = 16 
V, Va = Vs = V, Ve = Vf = 1 V, and magnetic field 5 = 
T. (b) The Coulomb blockade oscillations. The first 5 peaks, 
labeled as P1-P5, are shown. The data is taken under the 
same condition as that in (a). 



ture. Other authors proposed that the quantum dots are 
from Si nanoparticles in their point contact sihcon-on- 
insulator (SOI) devices. Our device is free of nanoparti- 
cles, thus different from them. Therefore, it is most hkely 
that in our system the negative interface trapped charges 
create potential barriers along the one-dimensional (ID) 
point contact channel and form a quantum dot.l^^I These 
trapped charges are always present in the Si/Si02, and 
the current industrial standard provides an interface trap 
density of approximately 1 x 10^° cm~^. In all of the four 
SETs we have characterized, the single electron tunneling 
characteristics are similar, although the details of which 
are sensitive to the specific interface charge distribution. 
There are changes due to thermal cycling, but the main 
features discussed here are stable and can be reproduced. 
Because the top gate threshold voltage is measured to 
be ~ 9 V for the SET, using a parallel capacitor model 
and that the quantum disc has a diameter of 60 nm, at 
Vgi = 16 V the number of electrons in the quantum dot 
is estimated to be at most ^ 10. As the source-drain dc 
bias V21 is fixed at V, the SET displays the Coulomb 
blockade oscillations, as shown in Fig. ^h). The fol- 
lowing discussion will be focused on the magnetic field 
dependence of the first five peaks labeled as P1-P5 in 
Fig.ib). 

We first focus on the evolution of the ground state and 
the excited states near peak 4 in magnetic field. The mea- 
sured excitation spectrum is shown in Fig.mb).^ In Fig. 
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FIG. 4: (Color online) (a) Differential conductance G21 as a 
function of V21 and the side gate voltage Vc,d , near peak 4, 
at the magnetic fields of 3.3, 4.8, 6.2, and 8.8 T, respectively; 
Vgi = 16 V, Va = Vb = V, Vb = Vf = 1 V. (b) The 
corresponding excitation stripe taken at V21 = —4.7 mV with 
magnetic field B between 1 and 9 T. (c) Peak positions as a 
function of B are extracted from the raw data in (b). The 
peak positions are in chemical potential, E = eaVc,D, with 
an arbitrary offset. The straight lines are from linear fitting 
of the data. The ground states Eq are only fitted from 1 to 
3 T. Their slopes are labeled in the graph, (d) Schematic 
showing the evolution of single-particle energy levels Eq and 
El driven by the Zeeman effect, where Eq is the ground state 
with a spin-up electron, and Ei is the first excited state with 
a spin-down electron, depicted in solid lines. In a two-electron 
system, the B-evolution of the singlet and the triplet states 
follow Eo{B) and Ei{B), respectively. Therefore we also label 
the singlet state as Eo and the triplet state as Ei . The cross- 
ing between Eg and Ei, the singlet-triplet transition, occurs 
at Bt. 



|3ja), there is an excited state near peak 4 {Vc,d is about 
—2.2 V). Figure [4]ja) shows that this excited state moves 
toward the ground state as the magnetic field increases. 
When the source-drain dc bias V21 is fixed at —4.7 mV, 
and the magnetic field and Vc,d are swept, the excita- 
tion stripe near peak 4 is shown in Fig. W[h). With the 
increase of magnetic field (0 < B < 6 Tj, the distance 
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between the first excited state Ei and the ground state 
Eq decreases. As i? > 6 T, the state Ei becomes the new 
ground state. The maxima in Fig. |4 b) are fitted with 
Gaussians (as functions of Vc,d), and Vc,d can be con- 
verted into chemical potential, using E = eaVc,D, where 
a — 0.09, as defined in Eq. (jlj.^-^ After an arbitrary off- 
set, the resulting peak positions (in chemical potential) 
are shown in Fig. |4](c). The absolute energy position 
of a peak in the excitation spectrum is determined by 
practical experimental parameters; we therefore focus on 
the difference between the most important features, that 
is, the energy difference between the "ground state" Eq 
and the "first excited state" Ei. To first order, both 
Eq and Ei show apparent linear dependence on mag- 
netic field B. We therefore fit Eq{B) (in the range 1 T 
< B < 3 T) and Ei{B) (in the range 1 T < B < 9 
T) by linear lines. Straight lines for the guidance of 
the eye are shown in Fig. gc). The energy difference, 
Ei{B) - Eo{B), (= [-0.137- (-0.018)](meV/T)-B), is 
—0.119 meV/T. The Zeeman splitting in bulk silicon is 
0.116 meV/T (= qhbB, with g = 2, fj,B is the Bohr mag- 
neton.). It is clear that the energy difference is dominated 
by the Zeeman effect. Figure |4](d) illustrates the Zeeman 
effect and the ST transition, where we only consider the 
two electrons in the outmost shell. The lower lying elec- 
trons, if any, are ignored, for they do not influence the 
spin dynamics. For a one-electron system, the ground 
state is a spin-down state in a magnetic field. The Zee- 
man splitting has the linear magnetic field dependence, 
±gfiBB/2. In Fig. ^d), the spin of the second electron 
is depicted by solid black arrows. When this second elec- 
tron is added to the quantum dot, it would interact with 
the electron already in the dot and form either singlet 
or triplet states. As shown in Fig. |4]^d), quantum se- 
lection rule dictates the spin of incoming electrons, that 
is, a spin-up state for the singlet state Eq, and a spin- 
down state for the triplet state Ei. At low magnetic 
fields, < B < Bt, the singlet state has a lower en- 
ergy than that of the triplet state. While at a higher 
magnetic field {B > Bt), the triplet state Ei becomes 
the ground state. The spin configuration of the ground 
states of the two-electron system is therefore controlled 
by the magnetic field. Were the Zeeman splitting the 
only effect to be considered, the magnetic field induced 
ST transition should occur at about 14 T. This estimate 
is based on the data shown in Fig. |4|c), and that there 
is ~ 1.6 meV energy difference between Ei and Eq at 
zero field. When the size of the electron wave functions 
shrinks for increasing magnetic field B, the interdepen- 
dence of Coulomb interaction and single-particle states 
becomes important. Here wc only consider the first- 
order corrections due to the electron-electron interaction 
in this two-electron system. The chemical potential is 
^is{B) = E'o+ gf^BB/2 + Coo{B) = Eo{B) + ACoo(B) 
for the singlet state, where Coq{B) is the direct Coulomb 
interaction when both electrons are in the gound state, 
E'o + Coo(O) = Eo{0), and ACoo(S) = Coo{B) - Coo(O). 
For the triplet state, the chemical potential is firiB) = 



E[ - gfiBB/2 + Coi{B)-\Koi{B)\ = E,{B) + ACoi {B) - 
A\Koi{B)\, where Coi(B) and Koi{B) are the direct 
Coulomb interaction and the exchange interaction re- 
spectively when one electron occupies the ground state 
and the other occupies the first excited state, E[ + 
Coi(O) - |i^oi(0)| = ^i(O), ACoi(B) = Coi{B) - Coi(O) 
and A\Kni{B)\ = \Koi{B)\ - |Xoi(0)|.'2J Both the direct 
Coulomb interactions and the exchange interaction in- 
crease with increasing B, since the size of the wave func- 
tions decreases. For the triplet state, we hypothesize that 
ACoi(i?) is largely compensated by A|iiroi(^)li due to 
the apparent linear B-dependcnce of the triplet state; 
while the singlet chemical potential f^siB) increases with 
increasing B. So the electron-electron Coulomb interac- 
tion (including the direct Coulomb interaction and the 
exchange interaction) plays a significant role here and 
drives the singlet-triplet transition at a much lower mag- 
netic field between 4 and 6 T. There are some minor 
deviations from the linear i?-dependence of the triplet 
state, which could be from other effects in addition to 
this simplified picture. A detailed analysis requires the 
precise information of the electrostatic potential in the 
quantum dot, which is beyond the scope of this paper. 
The conductance of the singlet state drops when B > 4 
T, and it becomes invisible after the ST transition. This 
is due to the transport blockade effect, which will be dis- 
cussed later. 

The magnetic field dependence of the first five 
Coulomb blockade peaks is shown in Fig.jsja). The peak 
position and the peak amplitude are extracted from this 
raw data by fitting each peak with Eq. ([T|. The shift 
of peak positions (in Vc,d) can be converted into the 
change in chemical potential using the formula:A£'(i?) = 
eaiVS'^o'^iB) - V^f^^^XO)]. In the stability chart [Fig. 
|3ja)], from the slopes of Coulomb diamonds, a is about 
0.15, 0.088, 0.090, 0.090 and 0.061 for peaks 1 to 5, 
respectively.!^ These changes in chemical potential refiect 
the evolution of the ground state energies. Figure ^h) 
shows the peak spacing (in chemical potential) between 
successive peaks as a function of magnetic field. Using 
the energy difference in our analysis not only minimizes 
the uncertainty due to charge fiuctuation and the long 
term drift in analog electronics, but also helps in iden- 
tifying the role of spin in the addition energy spectrum. 
Based on the Zeeman splitting, four straight lines with 
slopes of either g/isi? or —gfisB, with g—2, are plotted. 
Since the second peak only appears between 2 T and 8 
T, we compare the first 5 peaks in this field range. Each 
curve is arbitrarily offset for clarity. Because the chemical 
potential differences between successive peaks show the 
alternate slopes of giiB and —gfiB, it is clear that the 
shifts of the ground states show an spin-down spin-up 
filling pattern and are dominated by the Zeeman effect 
at low magnetic field (B < 4 T). Note that such odd- 
even-odd-even alternate filling pattern is different from 
that due to isolated single dopants in silicon.'^^ The fact 
that the data are well explained by the Zeeman splitting 
is consistent with our earlier finding that the magnetic 
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FIG. 5: (Color online) Evolution of Coulomb blockade peaks 
with the magnetic field B. (a) Drain current of the first 5 
peaks (P1-P5) with Vac = 0.1 mV, Vai = 16 V, Va = Vb = 
V, Ve = Vf = 1- Each trace is offset linearly with B. (b) 
Peak spacing between successive peaks and (c) peak ampli- 
tude as a function of B are extracted from the raw data in (a) 
fitted with Eq. |T|. In (b), the peak spacing is in chemical po- 
tential and is offset for clarity. The straight lines have slopes 
QI-IbB or —q^lbB, assuming the Zeeman splitting with g — 2. 
In (c), the arrows in the squares indicate the spin configura- 
tions for peak 3 (P3) and peak 4 (P4). (d) Schematic of the 
magnetic field induced spatial separation of spin-up and spin- 
down states. In a magnetic field, the conduction band minima 
of spin-down and spin-up electrons will split due to the Zee- 
man effect with — q^bB /2. Electrons are fully spin-down 
polarized in the source and the drain leads near the quantum 
dot. (e) Schematic showing the transport for a two-electron 
system with spin-down polarized leads. When B < Bt, the 
ground state is a singlet, which only allows a spin-up elec- 
tron to tunnel through, i.e., the spin blockade state. When 
B > Bt, the ground state becomes a triplet, which permits a 
spin-down electron to tunneling through. The spin blockade 
is thus lifted. 



field dependence of energy levels (both the ground state 
and the excited state) is primarily driven by the Zee- 
man effect. For the data reported in this work, the mag- 
netic field is applied perpendicular to the sample surface, 
therefore, the orbital effect might be expected. However, 
all of the observed magnetic field dependence can be well 
explained by the Zeeman effect with a g-factor of ~ 2, 
leaving orbital effect not a significant factor in our sys- 
tem. This is caused by the specific shape of our quantum 
dot. When we squeeze side gates C and D, the electron 
wave function is elongated along the source-drain direc- 
tion, the orbital effect becomes negligible. This leaves 
the B dependence of the peak positions dominated by 
spin effects.'^ 

The SET is fabricated on a Si (100) wafer, so the valley 
splitting from the two-fold valley degeneracy is expected. 
It is generally believed in a Si (100) inversion layer that 
the valley splitting ranges from ^ 1.5 meV, increases 
with increasing electric field or m agneti c field, and also 
depends on the Si/barrier interface) ^ '' l ^^ l But detailed re- 
lationships between them are not clear yet. The valley 
splitting in Si quantum dots has not been systematically 
investigated in experiments, primarily due to the lack of 
the excitation spectrum data. Some speculate that it is 
0.35 ~ 0.46 meV.^® Others indicate that it is o f order 
of a few meV and greater even at i? = T.^^SEHl jn our 
SET, the spin-down spin-up filling pattern requires that 
the valley degeneracy is lifted, and the excitation spec- 
trum suggests that the valley splitting should be larger 
than 1.6 meV. This is probably due to the strong lateral 
confinement in our quantum dot. Further investigation 
is necessary to reach an affirmative conclusion. 

Finally, we turn our attention to the tunneling peak 
amplitudes in Fig.|5ja). In a lateral quantum dot device, 
an electron tunnels into and out of the quantum dot from 
the two-dimensional electron gas (2DEG) leads. As illus- 
trated in Fig. [5jd), the electron concentration will vary 
from the bulk value to zero near the tunneling barriers; 
the Zeeman splitting effect will cause the separation of 
the conduction band minima of the spin-up and spin- 
down electrons in a magnetic field, which will result in the 
different populations of the spin-up and spin-down sub- 
bands. In GaAs, spin polarized injection is due to edge 
states in a magnetic field perpendicular to the 2DEG.'^ 
This edge-state picture is probably not valid in our sil- 
icon devices, because the Shubnikov-de Hass oscillation 
minima are not zero up to 9 T in the test hall bar de- 
vice. However, magnetic field can fully spin polarize a 
2DEG with = 1-5 x 10" cm'^ for B = 5 T solely 
with spin effects in a high mobility Si MOSFET.22 23 
the test hall bar device, the Shubnikov-de Hass oscilla- 
tions show that the 2DEG is spin polarized at least at 
B = 9 T for n2d = 7.12 x 10" cm^^^ -^^ jg reason- 
able to expect that electrons will be polarized into the 
spin-down state in the 2DEG leads near the tunneling 
barriers above some magnetic field, as illustrated in Fig. 
|5jd). This effectively makes the spin-up electrons tunnel 
through a thicker barrier, so the injection is dominated 
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by spin-down electrons. 

In theory, peak amplitude is proportional to the tun- 
neling probability and depends exponentially on how 
much wave functions in the dot and in the contacts over- 
lap with each other. There are two different mecha- 
nisms governing the amplitude modulation of Coulomb 
blockade peaks in a magnetic field. One is the spatial 
mechanism when an electron tunnels into different spa- 
tial wave functions in a quantum dot, not related to spin; 
the other is the spin-blockade mechanism due to spin- 
polarized injection. The lower edge of the excitation 
stripe in Fig. ^h) reflects the evolution of peak 4 as a 
function of magnetic field in Fig. [sj^a). The conductance 
of the singlet state drops dramatically when B > 4 T, 
while the conductance of the triplet state doesn't change 
much, as shown in Fig. |4|^b) and Fig. [5]^ a). This dif- 
ference can be explained by both the spatial effect and 
the spin-blockade effect, since the singlet state requires a 
spin-up electron which occupies a smaller wave function 
at the center, while the triplet state requires a spin-down 
electron which occupies a larger wave function. However, 
the electron states corresponding to peak 3 (a spin-down 
single electron state) and peak 4 (a singlet state) should 
have similar (single-particle) spatial wave functions; their 
peak amplitude should have similar dependence on the 
magnetic field if there is no spin blockade effect. In Fig. 
[sj^c), the amplitude of peak 4 decreases much faster than 
peak 3 when the magnetic field is larger than 4 T. A rea- 
sonable explanation is the spin polarized injection from 
the leads. ■^^ For 0.1 mV small ac excitation between the 
source and the drain, only ground states can lie in the 
transport window. For a singlet state, only a spin-up 
electron can tunnel into and out of the quantum dot, 
the amplitude of peak 4 decreases dramatically owing to 
this spin polarized injection, and as shown in Fig. [sj^c), 
it occurs at i? > 5 T. This forms the spin blockade. Af- 
ter B > 6 T, the amplitude of peak 4 increases with the 
magnetic field, because the ground state corresponding to 
peak 4 changes from a singlet state to a triplet state. For 
the triplet state, the incoming spin-down electron experi- 



ences lower and thinner barriers, thus the higher conduc- 
tance. Since the amplitude of Coulomb blockade peaks is 
determined by the two-electron spin configuration in the 
dot with spin polarized leads, it can be used to distin- 
guish a singlet state from a triplet one. Monitoring the 
amplitude is thus a form of spin blockade spectroscopy.!^ 



IV. CONCLUSIONS 

We have investigated electron spin in a nanometer- 
scale quantum dot in silicon. Using the unique bilayer- 
gated structure, narrow point contact channels are elec- 
trostatically defined. Along the narrow channel, a quan- 
tum dot is formed by potential barriers, probably due 
to negative charges near the Si/Si02 interface. Tunnel- 
ing spectroscopy clearly demonstrates the single electron 
tunneling characteristics. There are two important spin- 
related features. First, applying excitation spectroscopy 
on a dot populated by even number of electrons, we have 
directly observed the ST transition. Second, based on the 
observed ST transition, we further use amplitude spec- 
troscopy, i.e., tracing the magnetic field dependence of 
Coulomb blockade peaks, to identify the spin blockade 
phenomena. With a spin polarized injection, we demon- 
strate that the amplitude spectroscopy can be used to 
detect a single spin. 
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